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ABSTRACT 


A method is derived for finding the linear response and 
loading transfer functions for the lateral aerodynamic case 
of airship flight through atmospheric turbulence. The func- 
tions obtained are in a form that can be applied to the 
various spectral analysis methods used to predict surviva- 
bility currently employed by designers. A numerical example 
using the USS AKRON (ZR-4) is presented. The results show 
that peak motion response and loading occur when the en- 
countered spectral component has a wavelength equal to the 
airship length, and that simple feedback of heading angle 
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NOMENCLATURE 


hull cross-sectional area 


total buoyancy force of the airship, 
o*g*volume 


hull bending moment about the vertical axis 
distributed along the hull longitudinal 
axis 


longitudinal characteristic length of the 
airship 


mean chord of fin 


matrix of coefficients from equation 51 


nondimensional aerodynamic force in the y- 
2 


ai rseraon,  C....=2¥/o0uU. Ss 
y fe) 
nondimensional aerodynamic rolling and yawing 


moments respectively, (C,, Cc.) = 2(L,N) /pU,°Se 


nondimensional aerodynamic lift 
Cy for the fins alone, no hull interference 
nondimensional shear force Cine 28 (1) /0U,*s 


nondimensional bending and twisting moment 


A ey 


respectively, (C 
nondimensional time derivative, (c/2U,)a( )/dat 
gravitational acceleration 


turbulence forcing functions (equation 59) 


body-fixed coordinate measured normal to the 
hull centerline (positive up) 


cm 


h location of the vehicle's mass center 


h location of the empennage-assembly's mass 
center 


Sears' function corrected for finite aspect 
ratios by Filotas [Ref. 16] 


moments of inertia about the x- and Zz-axes 
product of inertia w.r.t. x~ and z-axes 


nondimensional moments and product of inertia 
; a —3 

1g. BI. /pse 

imaginary operator, i = y-1 


hull potential cross-flow factor from Jones 
and DeLaurier [Ref. 6] 


axial and traverse apparent-mass coefficients 


control gain of fin normal force to vehicle 
azimuth angle 


axially-aligned body-fixed coordinate origi- 
nating at the nose 


axial location of the buoyancy center 


axial location of the hull-fin intersection 
point 


axial location of the fin's aerodynamic center 


axial mass-center location of the entire 
vehicle . 


turbulent scale length 
rolling moment 


mass of the entire vehicle, including internal 
air and gas 


mass of the empennage assembly 


yawing moment 


10 


~ 


2Pi Hees ef! So nolzevo!l aA 


— 4 ~~... Dette 


i{TaPeTT Lab ales 


BG t oO @4eg losspor 
-OmM ACR Ssae 


i i208 sixes 
- % +Cog 
.- o' 2200G, Ieixs 


7 s7O1 3=7/102-?6en Jalen 
wu cistiev 


spné!l eff. Jneludse4 


eum pediflos 


Sige” Slee — nes a 


<> 


=> 


GREEK SYMBOLS 


O 
1e) 


B 


vehicle angular velocity about the x-axis 


nondimensional value of p, p = (c/2U,)p 


nondimensional maximum value of 5 
vehicle angular velocity about the z-axis 
nondimensional value of r, f= (c/2U,)xr 
nondimensional maximum value of f 
auto-correlation function 

: : aJi3 
reference area of airship, (volume) 


hull shear force, normal to the centerline 


stabilizer reference area (planform area) 


time 
twisting moment, about the hull axis 
reference flight speed 


perturbation velocity of the vehicle's mass 
center in the y-direction 


nondimensional value of v, % = v/U, 
nondimensional maximum magnitude of ¥ 


horizontal velocity of the atmospheric turbu- 
lence 


body-fixed wind-aligned stability axes (x 
positive forward, y positive right, z positive 
down ) 


axes fixed in inertial space 


force in the y-direction 


reference aerodynamic relative angle of attack 


aerodynamic sideslip angle 
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Y horizontal nondimensional velocity of the 
spectral component 


i maximum value of y¥ 


stabilizer efficiency factor, from Jones and 


g DeLaurier [Ref. 6] 
u nondimensional mass, u = 2m/oSc 
E axial coordinate, measured from the nose 
fe) atmospheric density 
fo} turbulence intensity 
G nondimensional stability root 
> roll angle 
0) maximum value of 6 
meso (hz) power-spectral function for turbulence com- 

JJ ponent v 

g 

v azimuth angle 
y maximum value of wv 
W spectral component frequency 
Q wave number 
SUBSCRIPTS 
( is aerodynamic force or moment terms 
(dp buoyancy terms 
ar control terms : 
( ene mass center terms 
( ) emp term for entire empennage assembly 
( Vg atmospheric-turbulence term 
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hull term 


inertial term 


reference equilibrium value 


derivative 


derivative 


derivative 


derivative 


fin term 


thruster-rotor term 


derivative 
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Vv 


derivative w.r.t. Dv 


derivative w.r.t. 
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derivative w.r.t. D8 


nondimensional term 


derivative w.r.t. time 
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I. INTRODUCTION 


When a close look is taken of the history of airship 
flight, it becomes evident that turbulence is a prime cause 
for concern in the design of lighter-than-air (LTA) craft. 
The spectacular crashes of the airships Shennandoah, Akron, 
and Macon are perhaps the most obvious reminders of the power- 
ful effect the wind can have on these fragile vehicles. It 
is imperative that proper consideration be given to the gust 
response of an airship, both dynamically and structurally. 

At the time the great rigids were built, designers had 
only a cursory knowledge of turbulence and how to deal with 
it. Burgess [Ref. 1], for example, dedicates only one page 
to the subject, and that is mainly a warning that gusts en- 
countered in flight can place larger loads on an airship 
than any maneuver of which it is capable. He suggests using 
the standard gust analysis technique of the day--a "fixed- 
in-space" flight through a ramp shaped gust. This was the 
method used to evaluate all aircraft at the time, and pro- 
vided some measure of confidence that the structure would 
withstand the stresses of flight. Its principal weakness 
was that it took into account only the direct gust loads. 

Since that time, the methods used to analyze response to 
turbulence have been greatly improved. Over time, the ramp 


shaped gust gave way to other gust shapes, finally settling 
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on the (l-cosine) gust as standard [Ref. 2]. By the mid- 
fifties enough data had been gathered on the nature of the 
turbulent wind to derive some basic statistical relations. 
This, in turn, allowed the development of the power spectral 
method for aircraft load analysis (see, for example, Press 
and Meadows [Ref. 3]). The spectral method has been applied 
to many types of aircraft over the last twenty years. At 
present, both methods are used in analysis of flight struc- 
tures, the one giving the more conservative structure de- 
termining the final design. 

While the development of turbulence modeling techniques 
progressed, applications to airship technology were slower 
in coming. Because the success of heavier-than-air (HTA) 
craft made the slower airship economically less attractive, 
LTA research lagged. In the period between World War II and 
the Arab Oil Embargo of 1973, the only significant contribu- 
tion to the study of airships in turbulence was by Calligeros 
and McDavitt [Ref. 4]. This paper presented a method of 
analysis that allowed a stable airship to respond dynamically 
to both sinusoidal and (l-cosine) gusts. Thus, the inertial 
and aerodynamic reaction forces of gust encounter could be 
included in the model. Further, by using the sinusoidal 
representation for gust shape, it was possible to apply spec- 
tral methods to the analysis. 

The recent interest in LTA brought about by rising fuel 


costs has increased the research in all aspects of airship 
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flight. The spectre of the great rigids disintegrating in 
turbulence makes it imperative that the designer have ade- 
quate means to predict an airship's response to gust pene- 
tration. Current research is aimed at supplying that means. 
DeLaurier and Hui [Ref. 5] refined the technique of Calli- 
geros and McDavitt [Ref. 4], by introducing refinements to 
the aerodynamic cross-flow model [Ref. 6] and allowing for 
stability augmentation through pilot control input. The 
model allows statistical prediction of an airship's dynamic 
response and operational lifetime for various combinations 
of speed, altitude and control gain. 

DeLaurier and Hui's work (as well as most others dealing 
with this subject) concentrated only on the longitudinal 
aerodynamic case. This thesis proposes to apply their model 
to the lateral case, enabling the response to side-force to 
be calculated. Bending and shear in the horizontal plane, 
as well as twisting moment, can then be taken into account 


when predicting airship life expectancies. 
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II. THE TURBULENT WIND 


The motion of the atmosphere is very complex. Shearing 
stress between layers of different speeds and at the ground, 
thermals caused by solar heating, weather fronts, vortex 
shedding behind obstructions and aircraft, plus many other 
phenomena, all contribute to a velocity field that is most 
difficult to describe. Of the methods available, that chosen 
will depend on the purpose for which it is used. A goal in 
the design of any flight vehicle is safety, with performance 
adequate for the mission. This dictates using models giving 
reasonable estimates for the design parameters, and it may 
or may not be necessary to closely match the physical reality 
to do this. In any case, the designer must be familiar with 
the advantages and limitations of methods available in order 
to choose wisely. What follows is a brief review of some of 
the turbulence models currently in use and how they apply to 


airship analysis. 


A. THE DISCRETE GUST 

As mentioned in the introduction, the discrete gust model 
has been used to analyze aircraft for many years. It is 
especially good when response to the passage through a steady 
velocity gradient, such as a thermal, mountain updraft, or 


jet stream, is desired. The method has been improved 
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steadily until, as Etkin [Ref. 7] points out, it has attained 
a high degree of sophistication. 

Figure 1 shows the shape for a (l-cosine) gust, where Wa 
and qd are maximum gust velocity and distance along the 
flight path of this maximum. By varying these parameters, 
the gust severity can be controlled. The value for qa is 


prescribed by the Federal Aviation Administration (FAA) as 
2d. = 25¢ 


The size of Wa is dependent on airspeed and altitude, and is 
shown in figure 2 for three values of equivalent airspeed. 
The 25c wavelength was chosen because it historically couples 
with the short period pitch mode of a rigid aircraft to pro- 
duce the largest load factors. Calligeros and McDavitt [Ref. 
4] showed that, for airships, the maximum loads occur when 
the wavelength is equal to the airship length. 

The British dictate (ARB CAR CH D3-3) that the gust para- 
meters be chosen to produce the peak response with aircraft 
flexibility taken into account. In this way, the model is 


"tuned" to the aircraft, thus assuring a conservative design. 


B. RANDOM TURBULENCE 

Extensive measurements of the atmospheric velocity field 
have been made, and the techniques involved are well estab- 
lished and reliable [Ref. 2]. They show that the velocity 


vector is best characterized by a random function of space 
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and time that is, in general, non-homogeneous, non- 
stationary, and anisotropic. The exact function has not 
been developed due to its tremendous complexity. Until 
enough data is collected (if ever) to allow the precise 
formulation, certain simplifying assumptions must be made 
to enable flight vehicle analysis. 

One assumption that applies everywhere except in the 
planetary boundary layer (below about 1000 ft), is that the 
turbulence is 'homogeneous', that is, the statistics of the 
field do not vary through space. In the boundary layer, 
scale length and intensity are homogeneous in the horizontal 
plane, but not vertically. Another assumption made is that 
the turbulence is 'stationary', or statistically time con- 
stant. Over the time periods of interest to flight this 
approximation is quite adequate. Also, the turbulence is 
assumed to be isotropic (again, except in the boundary 
layer), making the statistics invariant with orientation. 

One last simplification used to model atmospheric tur- 
bulence is the 'frozen field' or Taylor's hypothesis [Ref. 
8]. The change in the velocity field perceived by an air- 
craft as it passes with speed Ue through the air is given by 


the. substantial derivative 
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Taylor's hypothesis states that for all but the smallest 
values of Uo: the second term dominates and the first may 
be ignored. The result is that the correlations and spectra 
reduce to three-dimensional functions of space only. Phy- 
sically, this means the velocity field is 'frozen' in time, 
and the changes are due only to displacement. 

In an effort to specify an acceptable lower limit on Us: 
for which the frozen field applies, Dobrolenskiy [Ref. 9] 
cites studies comparing records of turbulence spectra gath- 
ered by captive balloon and an aircraft flying nearby at the 
same time. Within the margin of error, the two are statis- 
tically quite similar, and he concludes the lower limit on 
ae is comparable to the convection velocity (for all prac- 
tical purposes, the mean wind speed). Etkin [Ref. 7] points 
out that the vehicle speed can be as low as one-third the 
wind speed for good results. Note that the only vehicles 
capable of less than this velocity are LTA and VTOL craft, 
and then only when they are convected downwind with the air- 
mass. In hover, or upwind flight, the hypothesis holds. 
For that small portion of the flight envelope in which it 
does not, the forces generated on an airship's structure are 
small and, therefore, do not present a problem. 

The techniques for dealing with isotropic, frozen turbu- 
lence are well known [Ref. 8]. For those not familiar with 


the mathematical background necessary to deal with the 
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subject in depth, .Chapters 2, 3, and 13.of Etkin [Ref. 10] 
will provide a good primer. 
1. Turbulence Organization 

Anyone who has seen leaves swirl on an autumn day, 
or watched someone blow smoke rings, has an intuitive under- 
standing that the motion of the atmosphere is not completely 
arbitrary. What happens at one location effects conditions 
at another. Fuidynamicists characterize this interdependence 
by using expressions relating the stress and strain in the 
fluid. The technique has been applied to turbulence [Ref. 
11] with some success in predicting the actual velocity field 
of a boundary layer type flow. 

For flight vehicle analysis, a more convenient method 
of specifying the velocity field is the spectral decomposi- 
tion of the three-dimensional homogeneous vector field [Ref. 
12]. Figure 3 shows an aircraft flying through a (two- 
dimensional) sinusoidal wave of shearing motion. The velo- 
city change from the mean is given (for the lateral gust 


component) by 
1(Q,x' + Qoy') 
dv (x',y') =e dc, (1) 


where Qs and 25 are the wave number components in the x' and 
y' directions respectively, and C5 is the complex amplitude 


of the lateral component. If the vehicle penetrates the 
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field with velocity Uo: the coordinates become x' = x + Uot, 
and y' = y in the body fixed system. Equation (1) then 


becomes 


dv, (x,y) = SUE N reek er de, (2) 
The air velocity over the vehicle is then periodic with wave- 
lengths (27/2) 9) and frequencies (2) 2U,/2). The total 
field is made up of the superposition of these spectral 
components, much as a Fourier series represents a random 
scalar. 
2. Probability Distribution and Spectra 

With the expression for a single spectral component 
available, the next difficulty is determining the probability 
distribution of the individual frequencies. The power spec- 
tral density of a time varying function, X(t), is defined 


(in terms of wave number) as 


c 
@(2) = lim = | ES Ds) ae (3) 
ANQ+0 ere 
T +00 


where $(2) is expressed in (€t/sec)*/(radians/ft), and T is 
the duration over which X(t) is measured. The value is usu- 
ally computed by taking the autocorrelation function R(T) 
[Ref. 10: Chapter 2], and performing a Fourier transforma- 


tion, thus 
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Figure 4 [Ref. 2] shows the spectra of three different 
meteorological conditions. While varying in detail, they 
all exhibit the same decreasing trend at higher frequency. 
The vertical dimension is a measure of the intensity of 

the turbulence at that particular frequency, and the square 
root of the area under the curve a measure of the overall 
Eas oust velocity [Ref. 13]. .As a point.of-practicality, 
only the area under the actual measured curve is included, 
because values of the spectrum at higher frequencies con- 
tribute little to the response of an aircraft. 

From analyzing large numbers of samples, it is ap- 
parent the probability distribution is non-Gaussian [Ref. 
14], with very high and very low values of intensity occur- 
ring more frequently than predicted by a normal distribution. 
However, the vast majority of values do fall on a Gaussian 
curve, so it is reasonable to use this assumption for most 
applications. This is very beneficial because, whereas 
Gaussian input to a linear system results in Gaussian out- 
put, response to non-Gaussian input is, in general, unknown. 
In the analysis of flight vehicles, linear system models are 
used extensively, making the assumption of normal distribu- 
tion most desirable. In order to account for the large 


gusts omitted by this type model, the (l-cosine) method can 
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be employed. This is the practice recommended by many cer- 
tifying agencies. 
Two Gaussian models in current use are the Dryden 


spectrum 


eee aot 
le Se a ay ae) = 
(2+ G2) 
and the von Karman spectrum 
eae Nees PET BOER Heh he 
veo £6 ge : (6) 
33 7 PAVE Ky 
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The first was developed to define turbulence spectra in wind 
tunnels. It is the simpler of the two, but not as accurate 
as the second. For that reason it is not used as much today 
as in the past, but is mentioned here due to its historical 
Significance and the large number of references to it in the 
literature. 

Today, the von Karman spectrum is used almost uni- 
versally. In equation (6), o is the rms turbulence intensity, 
and L is the "scale length"--a measure of the average eddy 
size encountered. Testing has shown that the model is a 
reasonable fit to all levels of rareulenee: Figure 5 is a 
plot of one set of experimental data along with the predicted 
turbulence spectra for severe storm conditions using various 


values for L. As can be seen, the agreement is quite good. 
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The values for o and L are variable and appear to be 
functions of altitude. In addition, two standard categories 
of intensity are defined--"storm" and "non-storm". Table I 
[Ref. 15] lists values for non-storm (b,) and storm (b,) 
intensities, as well as scale length, currently used by NASA 
for horizontal atmospheric flight. In this table, the values 
Py and D5 are the probabilities of encountering non-storm 
and storm turbulence, respectively, at the altitude speci- 
fied. Note that the values of L given for flight below 1000 
ft are representative values that are probably low, giving 
conservative (high) numbers of load exceedances per unit 
length of flight. These values are usable for structural 
analysis, but inappropriate for control system studies of 
flight simulation where the vertical inhomogeneity must be 
taken into account. 

Equation (6) is the expression for the transverse 
spectrum needed to analyze the longitudinal response to 
vertical gusts. When dealing with lateral aerodynamics, such 
as the analysis done in Chapter III, it is necessary to use 


the longitudinal spectrum 


o*f, Z 


$1, (2) = 12 ——________ (7) 
11 nm (1 + (1-3392%0)7]>/° 


where o and 2 are the same as in the transverse case. The 
term “longitudinal” refers to the variation in gust velocity 


parallel to the direction of the mean wind, whereas 
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"transverse" refers to the perpendicular direction (both 


vertically and horizontally), thus 


It should be obvious from this discussion that flight direc- 


tion relative to the mean wind becomes important when dealing 


with lateral aerodynamics. 
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III. METHOD OF ANALYSIS 


A. THE MODEL 

The aerodynamic model used in this analysis is the same 
as used by DeLaurier and Hui [Ref. 5], except that it is ap- 
plied for the lateral case. The assumptions are as follows: 


i) the vehicle is perfectly rigid, flying at a reference 
velocity Uy through a constant density o0 


ii) the motions are described by the lateral case only 
iii) the control provided by rudder deflection is linearly 
proportional to the yaw angle (y), that is, 
ACye = KoY 
This last assumption is the case for a helmsman using greater 
control the farther off heading he is perturbed, an assump- 
tion in keeping with operational practice. 

The turbulence model used is the von Karman spectra 
described in Chapter II. It is assumed that the turbulence 
is composed of horizontal gusts only--either Ug or Vg de- 
pending on airship heading relative to the mean wind direc- 
tion. (vy is shown in the analysis.) 

1. Forces from Turbulence Components 


From Jones and DeLaurier [Ref. 6], the normal force 


on a hull segment of length dé is (see figure 6) 


rae | 2 : 9) dA ar a en - 
Fn = x90, [K sin(28) cos| | dt dé + (Can sind sin|98| 2rdé] 
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K is the hull potential cross-flow factor [Ref. 6] 


where: 
8 is the angle between the hull centerline and US 


ris the radius of the hull segment. 


Differentiating with respect to § to obtain a perturbation 


equation gives 


= 2 K cos (28) cos| 


Ne ee ae in(2 
adF, = Zou, | 5 sin(28) sin|3 
7 
{ 
P| 


(Ca), (cosésin| 6| + sindcos|9@|)d6 2r dé 


By limiting the analysis to the lateral aerodynamic case 


only, (9) becomes (8), the sideslip angle, and F, becomes 
the undisturbed 


Yue the hull sideforce. Further, if (So), 
value of sideslip, is assumed to be zero--the usual case-- 


the above expression becomes 


oe = ‘2008 (28.) cos 52 
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The stabilizer forces are given by 


2 
(v_) 
Se <M S 
(Tos 7 PF $5 (Crg) 5 H(kg) ng —y Se 


where: H(k.) is the generalized Sears function as given 
by Filotas [Ref. 16] 


we 
k_ = = , the “reduced frequency" of the fin 


oO 


(vo) is the gust velocity at the fin mid-chord 


The propellers used to drive the airship produce a 
side force when acted by the turbulence [Ref. 17]. Each 
thruster contribution adds to the total force and moment 
produced, and can be described, for the jth thruster-rotor 


combination by the following: 


2 
Cie) 4c 
" fe) at 
(gla, = 9-3 Sty (Cyg)a, (11) 
12) 
alg, = gla, Bom > Bry) (12) 
(Ng)n, = (glans og - bry) (13) 


Equations (12) and (13) assume that the rotors are arranged 
symmetrically about the x-z plane, so that moments due to 


rotor offset in the y-direction cancel out. 
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i Aerodynamic Forces and Moments Due to Airship Motion 


dv () 


1 
dt 


ne 2 aA. wl) [ . 
(dv on = 0U K>==—dae& —~ + a Rae ues (14) 


5 ae U. +Uork 


where Koy ky are the horizontal and longitudinal apparent 


mass coefficients respectively and 


W{e) =v = r{(li., en ie 


For the fins we have 


(y é Fo ag Z yes ) Vs v 2 ee cr 2 oe eve | is 
s = s || aoe Tita ob Mencia eee Clac age 
where ee =v- r(lon - 1.) 
i a r(lon - 1.) 
(Ls = eS Wane sen Ie (16) 
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ret 5 Z fae | ) or 
Maia 25 Ss Cs |Cny}, Us Chen) 


The superscript ac indicates the quantity in parentheses is 
taken about the fin aerodynamic center. 


Thruster forces and moments are given by: 


Ee Oo 
=-p-5- Sr, (Cyg)n, 7 (18) 
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where 


(vip) 5 Sr pee sae, - p(h,, - Br) 


3. Inertial Reaction of Airship to Aerodynamic Forces 


and Moments 

The "forces" covered here are those arising as reac- 
tions to airship motion. They are, in general, the negative 
of the forces and moments that give rise to the indicated 
airship translational and angular velocities and accelera- 
tions so as to produce a state of dynamic equilibrium 
(Ff. -ina=-0). 


For the hull: 


ee ay (ey ta) (21) 
where y' =v+ =(1. - €) + ULF 

(aL), = eh + <a (22) 

(aN), = =ar, r+ di, op (23) 
For the empennage: 

(Ys = ie Ms aay 

a oe ee ee eee (25) 

(Ls = - (Ty) sP 3 (Io) 5F (26) 
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where divs ale. and di, are the moments and product of 
inertia respectively of the differential element under con- 
sideration, including all structure, air and gas contained 
in the airship. 


4, Bouyancy and Control Terms 


Referring to figure 7, the force due to bouyancy is 


given by: 
(Yon = -(gdm - pgAdé) sind cosa, 
d = roll angle 
i steady state angle of attack 


Differentiating to obtain a perturbation equation gives 


(d¥ i, = (ogAdgé - gdm) cosh. cosa, dd 
and letting 6 = 0 results in 


(dye = [ogAd&é - gdm] cosa, do (er) 
Control force is assumed to come from rudder deflec- 


tion, and acts through the fin aerodynamic center. 


U 
(Yi), =e —- Sk (28) 


5. Shear Force, Bending and Twisting Moment 
The hull's shear force at station (1) is obtained by 


summing the sideforce values from the nose, up to (1): 
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where (dY)) = (dyn + (dyn - (dv Jn = (dyn 
and (a) is the number of rotors forward of station (l). 


Likewise, the bending moment at (1), measured along 


the centerline, is 


Pe : : ; 
BM(1) = | (lon &) (dY) + | (AN Jn Ly 2, hada) [(¥g) y+ (Fee) 
0 0 i ; 
(30) 
Finally, the twisting moment at station (1) is 
1 i al 
TM(1) = -| Dom (5) (d¥), + | (dim), + (dhig) hb 
0 0 0 
a 
- fp) CY “ae ke § 1 
os t,) 0 gt, * | w) 0) (31) 


The term (dL is the torque contribution due to the center 


mg) h 
of gravity being offset from the central axis (see figure 7). 


It is calculated from 


(dlig) h = Aa dm > cosa, 


B. FLIGHT DYNAMICS 


1. Dynamic Stability 


The equations for Lateral Dynamic Stability are taken 


from DeLaurier et. al. [Ref. 18] and given below: 


ee + ajay - (B - mg) cosa.é = 2u(D8 + £) (32) 
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Ba eed A = 02,8 cose.go =I, DS - 1, ,De 
Dd = 6 + 2 tana, 
DW =f Seca, 
Cy aero = Cy 28 + Cy 3D + Cyt + Cy;br + Cy F + Cys 
ACn sero = Cn2s + CngDs + Cy_f + Cy pdr + CnpP + Cn:Dp 
Beiearo * Clg8 + C1 gD6)+ C12 + Ci,Dr + C148 + C1528 
ACy, = Koy 
ieee (ts = tom! ee (ts 7 tom! a 
Cc 3 Bf = Cc 
acy u [h. : Dom! cy a [h. : Dom! Kv 
In this analysis 6 = = = Vv, thus equations 
(39) become i 
ACs, + ACY. = (B - tig) cosa, ¢ = 2u(DO + f) 
ACn sero + ACh, - 8 cosa.? = T,,Dt - TDP 
eye. + AC), = 2.5 cosa,¢ = 1,06 - pf 
Dd =p +f tana, 
Dy = f seca, 
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(32) 


(33) 


(34) 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


(41) 


(42) 


to 


(43) 


(44) 


(45) © 


(46) 


(47) 


ACy sero = Cy3¥ + Cy 3Dv + Cy? + Cy,Df + Cy pf ~ CysDB (48) 
a ee + C13bVv + cit + Cy ,DE + C156 + C15PP (50) 
These equations are linear, and along with equations (40) 
through (42) can be written in matrix form as: 
v Poe 
r Dr 
(cl | 8 - [D] Dp = 
9 D¢ 
n~ | wn 
1 D 
vo i ae 
where: v is translation in the y-direction 
ris the rotation rate about the z-axis 
p is the rotation rate about the x-axis 
¢ is the roll angle 
y is the yaw angle 
The matrices [C] and [D] are given by 
[Cc] = Cy. (Cy -2u) Cy (mg-B) cosa. “i | 
ae (Lo-lom) 
Cneg Cny Cnp “X,B cosa, = Ke 
Ae +h | 
Cle Ci+ Clp “Z)B cosa, cm | 
c c 
0 tana, ih 0 0 | 
0 seca 0 0 0 
° _ 
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To solve this system of equations, assume a solution of the 


form 
e auPig thc Sraigg, S2ect 12" Re: 19% 
(52) 
seg oy y = verot 
which, when substituted into equation (51) becomes 
 ¥ 
R 
ie = ied] 1 /B_| <0 (53) 
? 
es 


where ¢ is the non-dimensional stability root. This is an 
eigenvalue problem. A computer solution is performed to find 
the eigenvalues (stability roots) and eigenvectors (model 
vectors) for the control gains considered. For the subse- 
quent load-response analysis, dynamically stable cases must 


be chosen. 
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2. Turbulence Forcing Functions 


As mentioned above, the model chosen for the turbu- 


lence is sinusoidal with Gaussian statistics, in particular, 


v_(é) cosa_) 
—S— = 1 expliut - ive —— (54) 
fe) fe) 
By using this expression in equation (9), integrating from 


pae-O' to € = ly (the hull/fin intersection) and adding the 
contributions of the fins and thruster-rotor combinations, 
the complete turbulence forcing functions for the airship 


can be found. That is, 


1h a 
Ys = | dg) +) gla, + (%g)g (55) 
0 Soa 


Likewise, the yawing moment about the nose is 


1 
h a 
- +| Ed(¥g), +) [lps (%g)oJ + 15 (¥Q) 
2 j=1 ; J 


Ngnose 


eee “Gnoge’ 7 tomy ay 


and rolling moment is given by 


ee ae: ee 4 (57) 


G.) (58) 
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and the non-dimensional equations can be expressed as 


Gy = Gy, ¥ Gy ae C1. Gn = Gn, (59) 

where y = Texp(iwt) = Texp (ikt) (60) 
cw 
ko 
2U, 


Gy, G1. and G are the turbulent forcing functions for the 
Yy Y ny 


vehicle. 
3. Motion Response Transfer Functions 
Using the functions of equation (59) as forcing 
functions on the right-hand Side of equation (53), and 


dividing through by y gives 


ioe ale] vie sent 

bf 

R/T ¢ 
ny 

P/T = | G1 (61) 
£ 

an ie 

vier | | 0 


— — a — 


Solution of this expression for specific reduced 
frequencies (k), or spectral wave numbers (2), and fixed 
stable control gains (ky), allows calculation of the ex- 
pressions necessary for the solution of distributed force 


loadings and moments, by means of the following expressions. 
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ae iF exp(ikt), 3 exp (ikt) 
fe) fe) 
== E = exp (ikt), > = - exp (ikt) (62) 
fe) 
Vv ® “~ 
) = w exp(ikt) 


C. LOAD RESPONSE TRANSFER FUNCTIONS 
Once the motion response of the airship is known, the 
load response transfer functions can be calculated. 
1. Turbulence Loading 
These may be obtained by substituting equation 54 
into equations (9) through (13) and dividing by y (equation 


(60)). This gives, for example, for equation (9): 


2 
(ay -) U 
beeen: «ee Aeon = exp (-in& cosa_)dé (63) 
r 2 : g 
where 2 = _ = 2k 
° Cc 


Table II gives the complete list of load response transfer 


functions. 


2. Motion Response Loading 


The aerodynamic-reaction loading may be obtained by 
sas 
Uy 
tions (14) through (20) with the corresponding motion- 


’ & |, etc., in equa- 
° 


replacing the motion variables, 


response transfer functions in equations (62). For example, 


equation (14) becomes 
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(64) 


The Inertial-Reaction, and Bouyancy loading transfer func- 
tions are similarly obtained from equations (21) through 


(27). For example, equation (21) becomes 


2 
(dY_) 2U a é A - 
SR A = = 9 R 7 ae = “2 OR anes } 
aD = SF + igu*|= - SF1g 78) | fam (65) 
i GS » a < 
and equation (27) becomes 
(dY_) ) 
ik + 


3. Shear Force, Bending and Twisting Moment Transfer 


Functions 
The shear force loading transfer function is ob- 


tained using equation (29): 


en. oe pe — (67) 
0 J=1 
where: ‘4 
ay _ (dy dn i (dyn " (dy , (dye) 5 
; r ik r r 
and 
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and finally the twisting-moment transfer function, from 


equation (31), is: 


1 1 1 
dL dL 
TM(1) _ -| Se aitatint (dL) . | (dling) h 
r cm h a : 
0 0. 'oe 
(ey OE 
- } hp,.| + + ——-1 (69) 
ss iad iia r 


An important check on the analytical model is ob- 
tained by ensuring the net Shear Force, and Bending and 


Twisting Moments equal zero. 


S (1) (Y_) le ae Bs aR oe 
h Fs Gosis eS... mS . cs . 9 (70) 
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The moments are evaluated to the empennage assembly's mass 


center, (lJ 5, so that 
. 
BM(1__) (v) Cr fe pep Se 
2S PSR ee ae il pl: Galen a8 as S| 
cm’ s s ; E ‘ m4 
(N_) (N_) 
+—5+—T5 = 0 (71) 
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(72) 
These equations may be non-dimensionalized as follows: 
= ae aad (73) 
EF ou, S 
Wee og em 1) 
= 5 = (74) 
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D. RESPONSE TO ATMOSPHERIC TURBULENCE 

Once the force and moment transfer function coefficients 
are known, the turbulence statistics can be applied to ob- 
tain estimates of airship lifetime and failure probability. 
When dealing with the lateral aerodynamic case, two different 
spectra must be considered--9, , and 255. This analysis will 
use the von Karman spectra as given by equations (7) and (6) 
respectively. As explained in Chapter II, the first is used 
when the flight direction is perpendicular to the mean wind, 
and the second when the direction is parallel. 

1. Root-Mean-Square Responses 

The root-mean-square response to turbulence of any 


system parameter can be obtained using its transfer function 
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multiplied by the spectrum (provided, of course, the spec- 
trum is Gaussian). For example, the rms shear force coef- 


ficient is 


SS td (76) 


Response to various conditions can be evaluated using appro- 
priate values for o and L in the equations for daa and 
choosing either 1 or 255 according to desired flight 
direction. 
2. Mission Analysis Method 

The mission analysis method is a technique for esti- 
mating a flight vehicle's probable lifetime. The method is 
based on the probability distribution of encountering tur- 
bulence in representative flight operations [Refs. 7 and 15]. 
It is assumed the total flight is a sum of Gaussian patches 
[Ref. 2]. The formula for calculating the number of ex- 


ceedences is 


‘ |x-x_ | i Se ] 
N(x) = )tN_|p,exp ref + p,exp ef || (77) 
of? 1 = 2 i 
b.A boA 
a 2 
where: x = maximum structural value of bending, moment 


coefficient at a given station, 
x = value of x in one-g level flight, 
N(x) = average number of exceedences of the indicated 


value of x per unit time 
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No = number of zero crossings of x per unit time 
A = UNE eel O14 Vers 
t = fraction of time in mission segment 
Pi, Po = probability values from Table I 
Di, b, = intensity levels from Table I 
also 
: 1/2 
N 1 (BE) ng | 
fo) 27 (BM) ns | 
where 
2.-)2/ 
(oU_“se | 
(BM) .¢ = | | 2), an (78) 
ae fae 
0 
Dae Be 
[ou “se e 
(BM), = |— an 28, do (79) 
2 j T 
0 


The probable lifetime is then [N(x)]7+ 


Shear force and twisting moment are analyzed in the 

same fashion using the appropriate transfer functions. 
3. Other Methods 

DeLaurier and Hui [Ref. 5] also include Failure- 
Probability analysis [Ref. 3] and "Mil-Spec Storm" analysis 
in their paper. There are others in existence, such as the 
Design Envelope Analysis [Reft. 7] that have been used for 
HTA flight, and are well documented. Because of this, they 


will not be included here, as techniques for using them are 
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the same once the transfer functions of response are known. 
The exact method used is up to the designer, based on his 


needs. 
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IV. NUMERICAL EXAMPLE 


In order to illustrate the lateral aerodynamic case 
developed in Chapter III, an example using the USS AKRON 
(ZR-4) is presented. The flight conditions chosen are: 


us = 123 ft/sec 


ALE =" 1000: .£t 
The velocity represents the maximum for the vehicle, and the 
altitude is typical of its operational range. In addition, 
a condition of neutral buoyancy (B-mg = 0.0) was selected. 
The geometry was taken from Freeman [Ref. 19] and the weight 
distribution from Woodward [Ref. 20]. With this information 
available, the inertial properties of the AKRON could be 
calculated using the method of Scholaert and DeLaurier [Ref. 
21] (see Appendix). The values obtained are shown in Table 
III. The Hull cross-flow and stabilizer efficiency factors 
are calculated using the method given in reference 6, and 


are: 


K 0.93225 


n 0.2600 


Ss 

The apparent-mass coefficients are from Munk [Ref. 22]. 
The stability derivatives are taken from DeLaurier and 
Schenck [Ref. 18], and shown in Table IV. With these, the 


control gain, and the inertial and geometrical properties, 
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equation 53 can be solved to find the stable roots. This 
was done in reference 18 and the results are shown in figure 
8. Mode 4 (indicated in the figure) is characterized by 
roll, yaw, and sideslip of equal magnitudes not unlike the 
dutch roll mode of a fixed wing aircraft. Mode 5 is one of 
equal and opposite 8 and w motions with small $6 perturba- 
tions. Mode 6 is a relatively high-frequency rolling motion, 
little affected by control gain. Modes 1, 2, and 3 refer to 
longitudinal aerodynamic modes [Ref. 18]. From this analy- 
sis, a control gain of 0.2 was found to provide the minimum 
stable condition. 

The forcing functions were next obtained using equations 
54 through 60, and are plotted in figure 9. The peaks in 
all the curves occur at a wave number of about .008. This 
corresponds to a wavelength equal to the length of the air- 
ship. lth is significantly smaller than the others, as 
expected, due to the smaller moment arm through which the 
sideforce works in producing roll. 

With the turbulence forcing functions, equation 61 was 
solved to obtain the motion response transfer functions. 
These are shown in figures 10 through 14. Control gains of 
0.2, 1.0, and 2.0 were used to illustrate the effect of its 
variation. The most significant feature of these responses 
is the peak at a wave number of .008. This corresponds, as 
in the forcing functions, to a condition where the spectral 


component wavelength exactly equals the airship length. This 
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is the result predicted by Calligeros and McDavitt [Ref. 4] 
for the longitudinal case. DeLaurier and Hui [Ref. 5] also 
obtained this result, but only for cases of higher control 
gain. For the lateral case, the control gain does not sig- 
nificantly change this peak, although, for yaw and yaw rate 
(J and y), and to a lesser extent roll (6), the response at 
lower wave numbers is reduced. 

Finally, the load response transfer functions were cal- 
culated. The results are shown in figures 15 through 17 for 
a wave number of .009, and a control gain of 0.2. Complete 
results are given in the appendix. For the most part, the 
results yield no surprises. The magnitudes of the load 
response follow the general trend of the combined motion 
responses, thus the peak loads occur at wave numbers near 
-008. The location along the axis of the peak load varies 
as the magnitude of the motions increases, shifting aft in 
the case of shear and twisting moment, and to the center for 
bending moment. Again, the lack of significant change with 


control gain is apparent. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The analysis presented is an extension of the work by 
DeLaurier and Hui [Ref. 5], and is subject to the same re- 
strictions. That is, it is limited to small perturbations 
in order to allow a linear analysis usable with power spec- 
tral methods, and its ability to make precise predictions of 
the loading when used with one of the methods that accounts 
for severe turbulence is questionable. Nonetheless, it is 
a valuable tool in understanding the response to an initial 
disturbance, and when employed as the aerodynamic input to 
the various statistical methods discussed earlier, it can 
yield important design and operational insight. 

The limited effectiveness of the simple control model 
employed, which is typical of someone cuing his response to 
a compass, was demonstrated. It is suggested that an exami- 
nation of the effects of roll control and yaw rate feedback 
be made. This would allow a decision as to the feasibility 
of using control to provide gust alleviation. As discussed 
by DeLaurier and Hui, control gain made a large difference 
in the expected lifetime of an ere when considering only 
longitudinal aerodynamics. Undoubtedly, for the lateral 
case, even the yaw control used in this analysis will con- 
tribute to increased survivability due to the reduction of 


loads at low wave numbers. Thus, the next step is to employ 
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the statistical methods to discover how much change is 
realized. 

The case of combined longitudinal and lateral motion 
needs to be studied. No aircraft ever built has ever 
managed to fly through turbulence that is strictly one 
dimensional, as is assumed for this analysis. Coupling the 
two cases would give a much better idea of the true action 
of an airship in turbulence. 

Finally, some means must be found to establish the vera- 
city of this model, as well as that for the longitudinal 
case. To the author's knowledge, no investigation of the 
actual response of an airship to conditions of known turbu- 
lence has ever been made. This is, of course, a difficult 
project, considering the limited number of airships cur- 
rently available if full scale tests are to be carried out. 
Wind tunnel investigations, made in the various oscillating 
flow tunnels available, would be helpful. Until some tests 
are done, however, this type of analysis must be considered 
only for its qualitative aspects as opposed to its quantita- 


tive predictions. 
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COMPUTER OUTPUT FOR THE NUMERICAL EXAMPLE 
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TABLE III 
GEOMETRICAL AND INERTIAL PROPERTIES 


OF THE USS AKRON (ZR-4) 


total volume = 7,382,400 £t° 


ce = 785.0 ft 


¢ = 37,914 £t* 


364.24 ft 
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Il 
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tt 


337286 Et 
P= "s03.01 ff 


mass = 17,039 slugs 


buoyancy = 548,642 lb 


I, = 38,685,100 slug-£t? 
I, = 471,799,000 slug-ft? 

2 
Bos = 102,129,000 slug-ft 
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TABLE IV 
STABILITY DERIVATIVES OF THE USS AKRON (ZR-4) 


neutral buoyancy, i =e £e/ sec, ALT =. 1000 £t 


Cy, = -0.7224 Cy, = -0.9863 
Cy, = 70.3418 Cys 0.0586 
Cy, -0.0648 Cy, = -0.0913 
Cn, = -0.1710 Cné 0.0293 
Cn. = -0.2352 Cyne = -0-0991 
Cny = -0.0150 Gt 0.0028 
C1, = -0.0322 Cy, = -0.0456 
C1, = -0.0153 C1, 0.0028 
C1, = -0.0066 Cys = -0.0042 
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Figure 1. The (1-Cosine) Gust Shape 
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Figure 2. Derived Gust Velocity for Gust Loads Formula 
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Figure 3. Elementary Spectral Components 
in Two Dimensions 
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Figure 4. Typical Power Spectra of Vertical Gust Velocity 
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Figure 5. Measured and Fitted von Karman Spectra 
of Vertical Gust Velocity from Severe Storm 
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Figure 7. Schematic of Buoyancy Forces and Moments 
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Figure 8. Lateral Root-Locus of the USS AKRON (ZR-4) 
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